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Abstract 

We investigate the electromagnetic cascade initiated by injection of very high energy 
photons in the early Universe and calculate the cascade spectrum by solving a set of 
Boltzmann equations numerically. In the calculation we take account of Compton 
scattering off the background electrons and pair creation off the background nucleons 
as well as photon-photon processes and inverse Compton scattering. We also apply our 
cascade spectrum to the big bang nucleosynthesis with photo-dissociation processes due 
to heavy unstable particles and obtain the constraint on their lifetime and abundance. 
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I. Introduction 



Electromagnetic cascade initiated by very high energy electrons and photons is very im- 
portant in a number of astrophysical situations (e.g., Zdziarski 1988 and references therein). 
For example, very high energy 7-rays emitted from Cyg X-3 propagate through the cosmic 
microwave background and cause the electromagnetic cascade. Electromagnetic cascade also 
expected in active galactic nuclei, compact sources and galactic halos. In the early Universe 
the high energy photons (and electron-positrons) might be produced by some exotic sources 
such as primordial black holes and unstable heavy particles. 

The high energy photons emitted from the above sources scatter off the background 
photon and produce electron-positron pairs if their energy is greater than mg/e-y {nie'. 
electron mass, e^: energy of the background photon). Then the electrons and positrons 
produce the high energy photons by inverse Compton scatterings. This cascade process 
continues until the energy of the photons becomes too low to produce electron-positron 
pairs. When relevant processes are electron-positron pair production by photon-photon 
interaction and inverse Compton scattering, the spectrum of the cascade photons are found 
to be oc €^^'^ (e.g.. Berczinskii 1990). The spectrum is modified by the scattering between 
the high energy photons and the background photons as noticed by Svensson and Zdziarski 
(1990), and the spectrum obtained by Svensson and Zdziarski (1990) is good enough for 
most of astrophysical applications. However, when one study the electromagnetic cascade 
in the early Universe (z ^ 10"^), one should take other processes into account. In the early 
Universe, the number densities of electrons and nucleons are so high that the Universe is 
opaque to photons by Compton scattering off the background electrons and pair creation 
by the background nucleons. These Compton scattering and pair creation are main energy- 
loss processes for the photons whose energy is too low to produce electron-positron pairs 
by scattering off the background photons. Therefore these processes which are neglected in 
the previous work play an important role to determine the low-energy region of the cascade 
spectrum. 

In cosmological applications, it is essential to know the precise form of the cascade spec- 
trum. One of the most important cosmological effects due to electromagnetic cascade is the 
photo-dissociation of light elements (D, "^He, "^He). It is well established that D and He are 
synthesised in the early Universe when the temperature is between IMeV and lOkeV. The 
theoretically predicted abundances of these light elements are in good agreement with the 
obserA'ed ones. Howei'er this success of the big bang nucleosynthesis (BBN) is destroyed if 
the high energy photons produced in the cascade destroy the light elements. By considering 
this effect, constraints are imposed on energy sources such as massive neutrinos (e.g. Kolb 
and Turner 1990 and references therein). When one places these constraints, it is necessary 
to use the precise cascade spectrum. 
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In this paper we study the electroiiiagiietic cascade m the early Universe and obtain the 
spectrum of the high energy photons and electrons including the interaction with background 
electrons and nucleons (i.e. Compton scattering and pair creation). We have already applied 
our spectrum to the specific problems (Kawasaki and Moroi 1994a, 1994b). Here we describe 
how the spectrum is obtained in detail and apply the spectrum to BBN with general unstable 
heavy particles. 

Recently, Protheroe, Stanev and Berezinsky (1994) studied the similar problem by Monte 
Carlo calculation. Our method to obtain the spectrum is different from theirs and may be 
much more simple since our method is just to find steady- state- solution to the Boltzmann 
equations numerically. We also present some simple formula which fit to our cascade spec- 
trum. 

In Section II the relevant processes in electromagnetic cascade are reviewed. In Section III 
the method to obtain the cascade spectrum is described in detail. The application to BBN 
and constraint on the general decaying particle are given in Section IV. Section V is devoted 
to conclusions. 

II. Electromagnetic Cascade Processes 

Once high energy photons are emitted into the thermal bath in the early universe, 
they induce electromagnetic cascade processes. At first, we briefiy review the elementary 
processes which determine the shape of the spectrum. 

• If the energy of the high energy photon is sufficiently large, the phot on- phot on pair 
creation process (7 -|- 730 ^ e+ -|- e^, where 730 represents the background photon) is 
the most dominant process since the cross section or the number density of the target is 
much larger than other processes. However, this process has threshold, i.e. the energy 
of the initial high energy photon has to be larger than fn^ 1^^, with e^, being the energy 
of the background photon. As we will see later, this process determines the shape of 
the spectrum of the high energy photon for ^ m1j22T. 

• Below the effectii'e threshold of the double photon pair creation, high energy photons 
lose their energy by the photon-photon scattering (7 -|- 730 ^ 7 + 7). But in the limit 
of 0, the total cross section for the photon-photon scattering is proportional to 
and this process loses its significance. 

• Finally, photons lose their energy by scattering off' charged particles in the thermal bath. 
The dominant processes are pair creation in the nuclei (7 -|- N ^ -\- e~ -\- N ^ where 
N represents the nuclei in the background), and Compton scattering off the thermal 
electron (7 -|- e^Q 7 -|- e^, where e^Q represents the electron and positron in the 
background). 
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Emitted high energy electrons and positrons lose their energy by the inverse Compton 



scattering off the background photon (e"^ + 7bg 



+ 7)- 



Taking these processes into account, the Boltzmann equations for the photon and the 
electron distribution functions fj and /g can be written as 
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where DP (PP, PC, CS, IC, S, and EXP) represents double photon pair creation (photon- 
photon scattering, pair creation in nuclei, Compton scattering, inverse Compton scattering, 
the contribution from the high energy photon source, and the effects of the expansion of the 
Universe). Full details arc given in Appendix. 

Our purpose is to solve the above Boltzmann equations (1) and (2), and to obtain the 
high energy photon and electron spectra. In eqs. (1) and (2), there are two time scales which 
govern the evolution of the spectra, ^. e. the one is that of the electromagnetic processes and 
the other is the Hubble time. Since the former is much shorter than the latter, it is difficult 
to determine the evolution of the spectra by the numerical integration of the full Boltzmann 
equations. In order to solve eqs.(l) and (2), we use another method. When the photon 
source is active, the time variation of the source is small in one Hubble time for most 
cases. Then, the photon source can be regarded as stationary and the stationary state is 
established by the balance between absorption and emission processes at each epoch. Based 
on this fact, we have obtained the photon and electron spectra from the stationary solutions 
to the Boltzmann equations. In fact, to derive a stationary solution is much more practical 
than the full integration of the Boltzmann equations (1) and (2). In the next section, we 
will explain the basic methods which we have used. 



III. Numerical Integration of the Boltzmann Equation 

In deriving high energy photon and electron spectra, we hcive to use numerical method 
since the Boltzmann equations (1) and (2) are very much complicated. In this section, we 
explain our approach to the Boltzmann equations. 

At first, we classify the right-hand side of ccis.(l) and (2) into the "outgoing" parts and 
"incoming" ones. Below, we ignore the expansion terms (5/^/5t)|EXP and (9/e/5t)|EXP iii 
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the Boltzinaiiii equatioiib (1) and (2) biiice the expansion rate of the Universe is much smaher 
than the scattering rates of electromagnetic processes. Then, eqs.(l) and (2) can be written 
as 



dt 
dt 



= -re(ee;T)/e(ee) + /e,IN(ee), (4) 



where "incoming" terms /^jn and fe^m can be obtained as the sums of the contributions 
from the source terms and the functionals of the distribution functions of the photon and 
the electron; 



5/7(^7) 



roo 



dt 

POO 

+ / d<K,,e(e7,<;^)/e«), (5) 



+ 

S "'^^ 



/ de;Ke,^(ee,e;;T)/^(e;) 



dt 

+ r<A:e,e(ee,<;T)/e(e:). (6) 



The explicit forms of Ta and Kab (^7 B — j^e) can be obtained from the full details of the 
Boltzmann equations given in Appendix. 

As mentioned in the previous section, the strategy we use here is to calculate stationary 
solutions to the Boltzmann equations (1) and (2), which obey 

dt dt ' ^ ^ 

In this case, the distribution functions and /e of the photon and the electron can be 
formally given as 

= (8) 

The important thing is that the "incoming" terms /-,,iN(e,,) and fciNi^e) depend only on 
/-y(e') and /e(e') with e' > e^ , Cg. Therefore, if the distribution functions /^(e-y) and fe{^e) 
with e^, > e (and the source terms) are known, we can obtain /^(e) and /e(e) from eqs.(8) 
and (9). By using this fact, we have solved the Boltzmann equations (1) and (2) with 
monochromatic high energy photon injection. As one can see, extensions to the cases with 
non-monochromatic sources or the high energy electron source are trivial. 
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In our numerical calculations, we have used the mesh points Ei (0 < i < A/'mesh + 1) with 
Eq — Cmin and Ejy^^^^^ — Co, where e^nm is some minimum energy we are concerning (which 
we take e^nm — IMeV) and is the energy of the primary injecting photons. The energy 
range between ey^m and cq is divided logarithmically, i.e. the i-th mesh point Ei is given by 

Ei^e^iJ^] . (10) 

^ ^min ' 

By using the fact that the incoming photons into the mesh point E^^^^^^ are supplied only 
by the photon source, incoming terms for i — Ny^esh is given by 

^ ^ mesn 

feMEN^^J ^ 0, (12) 

where n.^ is the production rate of monochromatic photons, and A, = (-Ej+i — £'j_i)/2. 
Combining eqs.(ll) and (12) with eqs.(8) and (9), we can obtain the distribution functions 
at the A'mesh-th mesh point /^{E^^^^J and /^{En^^^J. 

Next we determine the distribution functions for lower energy region. Essentially, "in- 
coming" terms for the i-th mesh point are derived from the distribution functions fj{Ej) 
and fe{Ej) with j > i. Discretizing eqs.(5) and (6), "incoming" terms are (approximately) 
given by 

f-fMEi) ^ E {K-r-riEi. E,)f,{E,) + K,,{Ei,E,)UE,)} , (13) 

i>i 

f,,m{Ei) ~ Y.^i{I<n{Ei.E,)f^{E,) + K,,{E,,Ej)f,{Ej)}, (14) 
i>i 

from which we can obtain f-y{Ei) and fe{Ei) by using eqs.(8) and (9)."*^ In a way explained 
above, we have calculated the photon and electron distribution functions fj{Ei) and fe{Ei) 
at each mash points from i — Ny^esh — 1 to i = in order. 

In Figs. 1-3, we have shown the photon spectra for the cases e^o — lOTeV, ITeV, lOOGeV 
and T = leV, lOOeV. As we will see, distribution function f-d^-f) with e^^mg/22T is 
very much small compared with that of the low energy region, since the double photon 
pair creation process is effective for such high energy region. As we mentioned before, the 
scattering rate of the double photon pair creation process is much larger than those for other 
cascade processes (if kinematically allowed), and hence the photon spectrum are extremely 
suppressed for the region ^ m1j22T. Below this effective threshold, the double photon pair 
creation process receives a Boltzmann suppression and other processes become important. 
In fact, the photon-photon scattering is the most dominant process just below the threshold. 



^In fact, eqs.(ll) - (14) receive corrections of tlie order of A, which is due to tlie chscretization of the 
energy range. In our numerical calculations, we have taken these corrections into account. 
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while the shape of the photon spectrum for the sufficiently low energy region is determined 
by the scattering processes off' the charged particles in the background. 

As one can see in Figs. 1-3, both the photon spectrum for sufficiently low energy region 
and that for the photon-photon scattering region obey power-law spectrum; /-.(e-y) oc . 



We fit the photon spectrum at the sufficiently low energy region (e^ ^ mj/80T) as 

(15) 

and that for the photon-photon scattering region (m'^/SOT ^ e-y ^ mg/22T) as 



with being the total amount of the energy injection from the high energy photon source; 

5/7(^7) 



— J d6j6j 



dt 



(17) 

s 

The amplitude of the spectrum is proportional to the mean free time of the photon. In 
low energy region the mean free time is determined by Compton scattering and depends 
on the background temperature as ~ T~'^ . For the photon-photon scattering region, the 
amplitude is proportional to ~ since the cross section depends on ~ T"^. For the cases 
T = leV, lOeV, lOOeV and cq = lOTeV, ITeV, lOOGeV, lOGeV. we have calculated the 
fitting parameters Piow, -Ppp, ^^low aiid iVpp, and the results are shown in Table 1? These 
variables slightly depend on the background temperature, while their dependence on the 
initial photon energy cq is insignificant. 

Before closing section, we comment on the case with high energy electron sources. By 
numerical calculations, we have checked that the spectra in the case with high energy electron 
injection are almost the same as those with photon injection if the background temperature 
and the total amount of the energy injection are fixed. Numerically, their differences are 
at most 0(10%). Therefore, the formulae given in eqs.(15) and (16) with Table 1 are well 
approximated ones if one redefines pm given in eq.(17) as 

5/7 (V 

s 



Pin — J de^ej — — — + J de^ee 



dt 



(18) 



IV BBN and Photo-Dissociation of Light Elements 

The injection of high energy photons might destroy the success of BBN since the photons 
produced through the electromagnetic cascade destroy the light elements (D, "^He, "^He). 



^Notice that the formulae (15) and (16) are relevant for the cases when the initial photon energy eg is 
much larger than the effective threshold of the double photon pair creation process. 
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When the energy of the high energy photon is relati^'ely low, i.e. 2MeV e^, 20MeV, 
the D, T and "^He are destroyed and their abundances decrease. On the other hand, if the 
photons have high energy enough to destroy "^He, it seems tliat sucli high energy photons only 
decrease the abundance of all light elements. However since D, T and "^He are produced by 
the photo-dissociation of "^He whose abundance is iimch higher than the other elements, their 
abundances can increase or decrease depending on the number density of the high energy 
photon. When the number density of the high energy photons with energy greater than ~ 
20MeV is extremely high, all light elements are destroyed. But as the photon density becomes 
lower, there is some range of the high energy photon density at which the overproduction 
of D, T and "^He becomes significant. And if the density is sufficiently low, the high energy 
photon does not affect the BBN at all. 

From various observations, the primordial abundances of light elements are estimated 
(Walker et al. 1991) as 



where /)4He and are the mass densities of "^He and baryon. The abundances of light ele- 
ments modified by the high energy photons must satisfy the observational constraints above. 
In order to make precise predictions for the abundances of light elements, we have modified 
Kawano's computer code (Kawano 1992) to include the photo-dissociation processes. 

As an example of the high energy photon source, we study the generic unstable particle X 
which decays into a photon and a daughter particle {X 7-I-X') (Lindley 1985, .Juszkiewicz, 
Silk and Stebbins 1985, Kawasaki, Terasawa and Sato 198C, Ellis, Nanopoulos and Sarker 
1985, Ellis et al. 1992), and calculate the abundances of light elements with taking the mass 
mx, lifetime tx and the abundance Yxo = nxo/ {C{'i)T'^ /t^"^) of the unstable particle as free 
parameters {nxo is the number density of the heavy particle before decay). In the BBN 
calculation, the baryon-to-photon ratio rjB is also a free parameter. However as shown in 
the previous paper (Kawasaki and Moroi 1994a), the baryon-to-photon ratio 775 is not an 
important parameter because the allowed value for 775 is almost the same as that in the 
standard case. Therefore we fix 775 = 3 X 10"^*^ in this paper. Assuming that the unstable 
particles are almost at rest when it decays and that the emitted high energy photons have 

^In our analysis, we have not considered the constraints from the abundance of ^Li, since the cross section 
data for ^Be is not available and hence we cannot predict the abundance of ^Li a part of which come from 




(19) 




(20) 




(21) 
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unique energy eg = mx/2, {dfj/dt)\s can be written as 




(22) 



Then, comparing the calculated abundance of the light element with (19)-(21), we can obtain 
the constraints on the mass, lifetime and abundance of the unstable particle. (For details, 
see Kawasaki and Moroi 1994a.) 

In Figs4 - 6 the constraints are shown for the case cq being (a) lOMeV, (b)lOGeV and 
(c) lOTeV in the eoYxo ^ Tx plane. Since the spectrum of the cascade photons depends only 
on the energy injection for large cq, the constraint for cq = lOGeV is the almost the same 
as that for eg = lOTeV. As is mentioned before, if the lifetime is short (tx ^ 10^ sec), the 
photon-photon process is effective below the threshold of ^He photo-dissociation and only D 
and "^He are destroyed. On the other hand, for long lifetime, the photon-photon process is 
not so effective that the cascade photons destroy the ^He and result in the overproduction of 
D and ^He. When the energy of the injected photon is less than about lOMeV, only D can 
be destroyed (Fig. 4). In any case, the obtained constraints are very stringent, in particular, 
when the cascade occurs at t ~ lO'^sec or later. 

V. Conclusions 

We study the electromagnetic cascade in the early Universe. The spectrum of the 
cascade photons is obtained by finding steady-state-solution to the Boltzmann equations 
which describe the time evolution of the cascade photons and electrons. As shown in Table 1, 
the resultant spectra both in low energy region and in the photon-photon scattering region 
fit to the power low spectra (15) and (16), and the spectrum index in the low energy is —(1.2 
- 1.5), which indicates that the cascade spectrum is slightly softer than the spectrum with 
neglecting Compton scattering off the background electrons. This softening of the spectrum 
may be due to the dependence of Compton scattering on the photon energy. For high e^, the 
cross section of Compton scattering is proportional to e^^. This means that the mean free 
time of the high energy photons becomes longer as their energy increases. Since the amplitude 
of the spectrum is roughly proportional to the mean free time, Compton scattering relatively 
increases the amplitude in the high energy region and softens the spectrum. 

We also apply our cascade spectrum to BBN calculation with unstable heavy particle 
and obtain the stringent constraints on the lifetime and the abundance of the heavy particle; 



eoYxo < 10-^ - lO-^GeV 
eoYxo < 10-i3QgV 
eoYxo< 10-9GeV 



for lOSec < rx < lOSec, 

for 10*^sec<rx, eo>20MeV, 

for 10*^ sec < ry , cq < 20MeV. 



(23) 
(24) 
(25) 
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Here we places the constraints on unstable heaA'v particles. However, if we read the lifetime 
as the epoch for the injection of the high energy photons, the above constraints apply to 
generic photon sources. 

Finally we briefly coininent on the recent related work. Protheroe, Stanev and Berezinsky 
(1994) calculated the cascade spectrum including Compton scattering and other processes 
by using complicated Monte Carlo simulations. Since they do not present any spectra, we 
cannot directly compare our spectra with theirs. However they also obtain constraints on 
heavy unstable particles, which are consistent with our results. 



Appendix ; Boltzmann Equation 

To make this paper selfcontained, we show the details of the Boltzmann equations which 
govern the electromagnetic cascade processes in this appendix. 

(I) DOUBLE PHOTON PAIR CREATION [ 7 + 7 ^ e+ + e" ] 

For high energy photon whose energy is larger than ~ mg/22T, double photon pair 
creation is the most dominant process. 

The total cross section aop for the double photon pair creation process is given by 



aop (/?) = -Trri [1 - P') I (3 - P') In ^ - 2/? (2 - /?^J | , (26) 

where — a/nie is the classical radius of electron and (3 is the electron (or positron) velocity 
in the center of mass frame. Using this formula, one can write down [dfj/dt)\Dp as 



dt 



1 1 f°° 1 - 



4m 



■ (27) 

/3=Vl-(4™l/.) 

The spectrum of final state electron and positron is obtained in Agaronyan, Atyan and 
Najapetyan (1983), and (9/e/9t)|DP is given by 



dt 



= ]^r!< r de,^ r de,%^G(e„ e„ e,), (28) 



DP ^ ^= 

where represents the distribution function of the background photon at temperature T, 

t 1 



fji^j) — 2 f- /rr\ 1 7 (29) 

TT^ exp(e^/i ) — 1 



and function G{€e,€.y,€.y) is given by 

, 1 ^7) 



77 ^7; 

4 



ee< ml{e^ + Q \ (cg + <) J e^^^ 2 
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with Cg = + — Cg. 

(II) PHOTON-PHOTON SCATTERING [7 + 7^7 + 7] 

If the photon energy is below the effective threshold of the double photon pair creation, 
photon-photon scattering process becomes significant. This process is analyzed in Svensson 
and Zdziarslii (1990) and for e' ^ O^ml/T), {df-Jdt)\pp is given by 



dt 



35584 2 2 -6 /'"^ , r / A 2 fi ^7 (^-\'\ 1- -Mr- ^ 

= 10125^^"^^"^^ k ''''^''^'''^'-'Y-7, + [7j I I ^^-S/tK) 

1 946 , , roc _ 

-Me.yrim:\' de^ej^ie,). (31) 



506257r' 



For a larger value of 7, wc cannot use this formula. But for high energy photons, photon- 
photon scattering is not signihcant because double photon pair creation determines the shape 
of the photon spectrum. Therefore, instead of using the exact formula, we have taken m^/T 
as a cutoff scale of (5/-y/5t)|pp, i.e. for e'^ < ni-i/T we have used eq.(31) and for e'^ > inijT 
we have taken 



dt 



0. (32) 



pp 

Notice that we have checked the cutoff dependence of spectra is negligible. 

(Ill) PAIR CREATION IN NUCLEI [7 + iV^e++e+iV] 

Scattering off the electric field around nucleon, high energy photon can produce electron 
positron pair if the photon energy is larger than 2me. Denoting total cross section of this 
process apc, {df-y I dt)\^Y> is given by 



dt 



= -nN(TYc{<^^,)f^Xe.X (33) 



NP 



where is the nucleon number density. For crpc, we have used the approximate formula 
derived by Maximon (1968). For near the threshold (e^ < 4me), the approximate formula 
is given by 

. M 2 ^2ti (k-2V f 1 23 11 . 29 4 ,A . ^ 
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where 



rrip 



2k -A 



P = 



k + 2 + 2V2k 



(35) 



and Z is the charge of nuclei. For large (e^ > 4me), apc are expanded in the parameter 
k~^, which is given by 



^2 2 

Z ar^ 



28 , , 218 

— In 2A; 

9 27 



+ 1^ (In 2^)3 - (In 2^)2 + " t) ^""^^ + ^^^^^ + 

-(f)'{^--«} 

} + 0(A;-'^) 



7r^_ 7' 
6 2 



2\*^ r 29 , , 77 

- <^ In 2k 

kj 12304 13824. 



(30) 



Differential cross section for this process dapc/ dEj^ is given in Berestetskii, Lifshitz and 
Pitaevskii (1971); 



dapc 2 2 ( P+P- 



4 p p pl+pi 

- - ZE+E_ 2 2 

3 p+pl 



^P+P- 



1 



P+P- 



+ 



ElEl + plpl - m'iE+E. 



m% f E+E^pX E+E_ - pi 

{ + 15 h {_ 5 



where 



^P+P- 

P± 
L 



P^ 



p'l 



E'i. — m?. 



± 



In 
In 



E+E_ +P+P- + rn'j 
E^E_ — p+p- + ml ' 
E± ^p± 
E±-p± 



(37) 

(38) 
(39) 

(40) 



with E- {E+) being the energy of electron (positron) in final state. By using this formula, 
{dfJdt)\^P is given by 



dt 



f~ dapc 
— rijy I de^ 



NP 



(41) 
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(IV) COMPTON SCATTERING [7 + 6-^7 + 6"] 

Compton scattering is one of the processes by which high energy photons lose their 
energy. Since the photo-dissociation of light elements occurs when the temperature drops 
below ~ O.lMeV, we can consider the thermal electrons to be almost at rest. Using the total 
and differential cross sections at the electron rest frame acs aiid dacs/d^e, one can derive 



dt 



dt 



cs 



cs 



- n. 



6?acs(4'^7) 



dcTcsje-i + me - 4, e^) 
de', 



(42) 



(43) 



where rie is the number density of background electron. Explicit form of acs aiid dacs/de^ 
are given as 



1 



X 



( 4 8 A , , ,1 



1 



— Tir 



2 ^^^g j 



2m 



x 2(l-|-a:)2 
/I 1 \ 1 



7 "7 "7 



e \ / 



with 



s — ml 2e- 



rrif 



(44) 
(45) 

(46) 



(V) INVERSE COMPTON SCATTERING [ e± + 7 ^ e± + 7 ] 

Contribution from the inverse Compton process is given by Jones (1968), and {df /dt)\ic 
is given by 



dfJe^) 



dt 

mo 



dt 



ic 



ic 



= 2'Krlml 



de, 



2/e(ee) 



°° , /e(ee) 



J-p -^7(^7) pi ^ r -p \ 



O 2 2/e(^e) 

-27rr,m,^^ 
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where function -F(e^, Cg, e^) is given by 



F{e^,e,,ej)\o>g>i = 2g In g + (1 + 2g)(l - g) + , /^^'^j,' , (1 - g), 



2(i-r,g) 



otherwise 



0, 



(47) 



(48) 

(49) 
(50) 



with 



4e e 
mf 



re(^^e ^7) 
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e^o = 10 TeV 



Temperature 


iVi„„GeV^ Fpp %pGeV^ 


1 eV 
10 eV 
100 eV 


-1.57 1.6x10^ -5.10 G.9 X 10-1^ 
-1.34 5.4 X 10^ -5.20 6.0 x 10"^^ 
-1.22 1.7 xlO^ -4.84 1.1 x lO^^^ 


e^o = 1 TeV 


Temperature 


Plow iViow GeV^ iVpp GeV^ 


1 eV 
10 eV 
100 eV 


-1.56 1.4 X 10^ -5.07 6.2 x 10"^^ 
-1.34 4.9 X 10® -5.17 5.5 x IQ-^^ 
-1.22 1.4 X 10^ -4.79 1.0 x 10"^^ 


e.o = 100 GcV 


Temperature 


/'low iViowGeV^ Ppp iVppGeV^ 


1 eV 
10 eV 
100 eV 


-1.56 1.4x10*^ -5.01 5.7x10-1*^ 
-1.33 4.7 x 10® -5.15 5.3 x 10"^® 
-1.22 1.3 X 10^ -4.74 1.1 x 10"^^ 


e^Q — 10 GeV 


Temperature 


Plow iViowGeV^ Ppp iVppGeV^ 


1 eV 
10 eV 
100 eV 


-1.33 4.5 X 10® -5.12 5.5 x 10"^® 
-1.22 1.3 X 10^ -4.77 9.6 x 10"^® 



Table 1: Piow, ^low, Ppp iVpp for the cases of T = 1 eV, 10 eV, 100 eV, and eg = 

lOTeV, ITeV, lOOGeV, lOGeV. 
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Figure 2: Same as Fig.l except for cq — ITeV. 
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Figure 4: Constraints on the lifetime and the abundance of the unstable particle from BBN 
for Co — lOMeV. Top horizontal axis shows the redshift corresponding to Tx- 
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Figure 5: Same as Fig. 4 except for cq — lOGeV. 
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Figure 6: Same as Fig. 4 except for cq — lOTeV. 
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